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The thermodynamic distribution function for exclusion statistics is derived. 
^ . Creation and annihilation operators for particles obeying such statistics are 

^ i discussed. A connection with anyons is pointed out. 

P^ I Exclusion statistics, which is a generalized version of Fermi- Dirac statistics, was intro- 

duced a few years back |^. While physical particles in the three dimensional world obey 
either Bose-Einstein or Fermi-Dirac statistics, this generalized statistics is expected to find 
CK . application in regions of condensed matter where the number of spatial dimensions is in 

^ ! effect reduced. There has been an enormous amount of interest in the subject in the past 

one year 0-^, no doubt because of the wide range of applications envisaged. 

It was mentioned in |]I| that no microscopic theory was introduced there: a phenomeno- 
1^ • logical framework for appropriately detailed theories was presented. However, the basic def- 

inition of exclusion statistics, whereby the exclusion of particles from single-particle states is 
k> ! in general less limited than in Fermi statistics, tempts one to follow the standard discussion 

H ' of fermionic thermodynamics with the minor changes that seem to be necessary. Statisti- 

cs ' ri?Tr 

cal thermodynamics has recently been generalized in this spirit [Q. The distributions thus 
obtained have been made use of in 0,|]. 

Although the starting point of for the number of ways of putting particles in a 
collection of single-particle states represents an interpolation between the forms which are 
known to be valid for Bose and Fermi statistics, it is not supported by a direct combinatorial 
argument, except in the limit of large collections of single-particle states f^. However very 
large collections or cells cannot be considered together unless the temperature is very high, 
so departures from the distribution found in this way can be expected. Furthermore, one 
runs into negative probabilities even for simple, rational values of the degree of exclusion 
(P). It therefore becomes necessary to derive the generalized distribution without appealing 
to large cells. That is the purpose of this letter. 
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As mentioned above, the distribution found in has been used in 0, where the algebra 
of creation and annihilation operators for oscillators obeying this distribution has been 
derived. In the same way we shall discuss creation and annihilation of oscillators obeying 
the distribution derived in this letter, albeit in a much more limited way than [Q. An 
interesting connection with anyons will emerge. 

In the original formulation |jl]] , a degree g of exclusion may be said to exist if n particles 
in a single-particle state reduce the dimensionality of the Hilbert space available to an 
additional particle by a fraction gn. We have omitted here the indices introduced in [^ to 
accommodate several species. This degree is clearly and 1 for the normal bosonic and 
fermionic cases respectively. In general g has to be rational, but it is simplest when - is 
a positive integer, m say. Then one may define the exclusion condition by stating that a 
given single-particle state may accommodate only a maximum of m particles. This is the 
definition which will be used hereafter. 

The approach of ||^ is to approximate the number of independent ways of putting n 
particles in d single-particle states by 

[d + (I - mn - m 

W = ^ ^ — (1) 

This expression is exact for the bosonic and fermionic limits. One may seek to replace it 
for the general case by the number obtained directly from the definition of exclusion given 
above. But there is a simpler procedure avoiding combinatorics. 

The grand canonical partition function of a system of particles with chemical potential 
/i at a temperature f3~^ is 

Z = 5]e~'3(^-^^), (2) 

where the sum is over different states of the system as a whole with energies E and numbers 
A^ of particles. If the energy can be assumed to be additive, one can rewrite this as 

^ = ^g-/3E"»(^»-/^), (3) 

where the grand sum is now over the occupation numbers rii of the different single-particle 
states with respective energies e^. A state of the system is supposed to be uniquely specified 
by the occupation numbers of the individual single-particle states. The grand sum factorizes, 

Z = nEe~^"'^''"''^ (4) 

each factor corresponding to one particular single-particle state. Upto this stage everything 
is standard. The sum over rii is an infinite geometric series for bosons because all nonneg- 
ative integral values are allowed. For fermions, the sum has only two terms, because the 
occupation number must be either zero or unity. In the case of particles obeying exclusion 
statistics, the allowed values of rij are 0, 1, 2, ... m. The finite geometric sum over Ui can 
again be found: 

p-P{rn+l)(ei-^) _ i 



This grand partition function can be used to calculate all thermodynamic quantities. In 
particular, it is of interest to write down the average number of particles: 

N = ^[(e^(^^-^) - 1)-^ - (m + i)(e/3(-+i)(^>-M) _ i)-i]. (6) 

i 

As usual, this is interpreted to mean that the average occupation number of the i-th single- 
particle state is 

n, = [(e^(^'-'^) - l)-i - (m + l){e^("^+^)(^^-^^) - \)-\ (J) 

This can equivalently be written as 

_ _ g{m-i)/3(£.-/.) _^ _^_ + (m - l)e^('^~^) + m 

This generalizes the Fermi form in a straightforward way. At zero temperature, all single- 
particle states with energies below // are fully occupied i.e., have an occupation number 
exactly equal to m and all other states are empty, as in the distribution function of [Q. 
However, at finite temperatures, the functions are different. Note that the right hand side 
of @ lies between m and and decreases monotonically as e^'-'^'"^-' increases. 

It is instructive to consider the simplest case of new statistics, viz., m = 2. Here the 
average occupation number is 

g/3(<:i-M) _|_ 2 
'^' " g2/3(e,-/.) _^g/3(e,-/.) _^ ^ ' ^^^ 

This equation can be inverted by defining 

X = e^^^-'^) (10) 

and solving a quadratic in x in terms of fi: 



\/l + 6n — 3n^ + 1 — n , . 

X = . (11) 

This root is always real and positive because n cannot exceed 2. 

The method of ^ can be used to find the algebra of creation and annihilation operators 
corresponding to the distribution (pT]): 



aa^ = — a^a. (12) 

Here the standard notation for oscillator operators is used, so that A^ is now the number 
operator for the oscillator. This relation looks complicated, but it is easy to see that it 
corresponds to the matrix representation 

] 1 0\ 

(13) 




This is different from that found in on the basis of the distribution of [Q]. The number 
operator, which has the matrix representation 




N= \0 1 I , (14) 

can be expressed as 

N = a^a- aa^ + 1. (15) 

Although the operator algebra (|T^) is satisfied, the matrices also satisfy the much simpler 
q-oscillator algebra 

aa^ — qo)a = q~ , (16) 

with q = exp(— y). Thus the particles satisfying exclusion statistics of degree 1/2 appear to 
be associated with the q-value of not ±i as might have been naively expected, but exp(— y). 
q-commutation enters the discussion in another way if one considers different states of 
the system. The allowed range of occupation numbers of a single-particle state have been 
taken to be unaffected by the occupation numbers of the other states. This corresponds 
to product Hilbert spaces and essentially independent creation/annihilation operators for 
different single-particle states. One can write 

where Qij is a new phase factor. The phase factor for i ^ j is in general a matter of 
convention, but it turns out that a very natural choice is 

qij=exp{—eij). (18) 

To see this, note that the c^- symmetrized wavefunction of a three-particle state can be con- 
structed out of three single-particle wavefunctions //^ (with ^4=1, 2, 3) 0] as 



oi [ji Jj Jk + Qij Jj ii Ik + Qik Ji JkJj + Qik Qjk JkJi Jj + Qij lik Jj JkJi + %j Qik ^jk JkJj Ji J- 



1 

r,% [J i J ] J k ' 'iij J J J I J k ' 'ijk J I J k J J ' "ilk 'ijk J k J i J j ' 'iij "iik J ] J k J i ' "iij 'itk "ijk J k J j 

(19) 



This {—q) — determinant is in general nonvanishing, but if the three wavefunctions / are 



identical, it vanishes with the choice ([T8|) because of the property 

1 + g + g' = 0. (20) 

The generalized exclusion principle with tti = 2 is then extended to states obtained by 
superposing the original states labelled by the indices i. 

We have considered the simplest case m = 2 in detail. For higher values of m, the 
expression for n contains powers of x higher than 2 and more complicated equations have 



to be solved to obtain x and hence the algebras. We shall not do that, but must point out 
that ( |T8|) generalizes to 



The {—q) — determinant form of the (m + l)-particle wavefunction vanishes again if all 
the single-particle wavefunctions are the same. This provides more direct evidence than 
that anyons with these q-values obey exclusion statistics. 

To conclude, we have derived the thermodynamic distribution function for analogues of 
fermions where a single-particle state can accommodate a limited number of particles. This 
is expected to find applications in condensed matter physics. 
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